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"9 Abstract 

^ ! A general formulation of noncommutative or quantum derivatives for oper- 
as ; 

£h ators in a Banach space is given on the basis of the Leibniz rule, irrespective 

, of their explicit representations such as the Gateaux derivative or commutators. 
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Os ; This yields a unified formulation of quantum analysis, namely the invariance of 
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quantum derivatives, which are expressed by multiple integrals of ordinary higher 
derivatives with hyperoperator variables. Multivariate quantum analysis is also 
formulated in the present unified scheme by introducing a partial inner derivation 
and a rearrangement formula. Operator Taylor expansion formulas are also given 



Q-i by introducing the two hyperoperators 5a->b = —^Sb and cIa^>b = ^(-(j- 1 ^)-^ 
> ! A 
with the inner derivation 5a '■ Q | — > [A, Q] = AQ — QA. Physically the present 



x 
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noncommutative derivatives express quantum fluctuations and responses. 
I. Introduction 

Recently noncommutative calculus has attracted the interest of many mathe- 
matians and physicists 1-15 . The present author 10-15 has introduced the quantum 
derivative df(A)/dA of the operator function f(A) in the Gateaux differential 1-3 

dm = ^A + MA)-J ( A) ^ dA 
h->o h aA 
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Here the quantum derivative df(A)/dA is a hyperoperator 10-15 , which maps an 
arbitrary operator dA to the differential df(A) in a Banach space. There is also 
an algebraic definition 8 ' 9 ' 12 ' 13 of the differential df(A) as 

df(A) = [HJ(A)] (1.2) 

for an auxiliary operator H in a Banach space. This differential depends on H. 
In particular, we have 

dA=[H,A}. (1.3) 
The property that d 2 A = requires the following condition 

[H,[H,A]] = [H,dA] = 0. (1.4) 

In the previous papers 10-13 , we have shown that the differential df(A)/dA 
defined in (1.1) is expressed by 

df(A) _ S m 

where 5 a denotes an inner derivation defined by 

6 A Q=[A,Q] = AQ-QA (1.6) 

for an arbitrary operator Q in a Banach space. The ratio of the two hyperoperators 
in (1.5) is well defined 10-13 when f(A) is a convergent operator power series. 

On the other hand, the derivative df(A)/dA defined through the commutator 
(1.2) is also expressed 9 ' 12 ' 13 ' 16 by Eq.(1.5). This is easily derived as follows. From 
Eq.(1.2), we have 



df(A) = 5 H f(A) = -5 f{A) H = -5 A l 5 A 5 m H 

= 5- A l 5 m {-5 A H) = S^SfwlHtA] = 6 ~jfdA, (1.7) 

using the commutativity of 5a and The meaning of the formal inverse 5 J 1 

in Eq.(1.7) will be discussed in the succeeding section. The above results suggest 



that the quantum derivative df(A)/dA defined in Eq.(1.8) is invariant for any 
choice of definitions of the differential df(A). One of the main purposes of the 
present paper is to make a unified formulation of quantum analysis and to prove 
the invariance of the quantum derivative df(A)/dA defined in 

df(A) = ^ • dA (1.8) 

for any differential df(A) satisfying the Leibniz rule 

d(fg) = (df)g + f(dg). (1.9) 

In Sec. II, some mathematical preparations are made on the formal inverse 
52 of the inner derivation 5a- In Sec. Ill we present a general formulation of 
quantum derivatives using the hyperoperators S^SdA and 5^_ 6 -i dA y A . Theorem 
I states the invariance of the differential df(A) for any choice of definitions of 
df(A). Theorem II gives the invariance of the derivative, df(A)/dA. Theorem III 
presents algebraic expressions of higher differentials {d n f(A)}. Theorem IV gives 
multiple integral representations of higher derivatives {d n f (A) / dA n } . Theorem 
V presents a general Taylor expansion formula of f(A + xB) in terms of higher 
derivatives {d n f(A)/dA n } for the noncommutative operators A and B. A shift- 
hyperoperator Sa(B) : f(A) \— > f(A+B) is also formulated. A general formulation 
of multivariate quantum analysis is given in Sec. IV, by introducing a partial inner 
derivation and a rearrangement formula. In Sec.V, an auxiliary operator method 
is briefly discussed, and it is extended to multivariate operator functions. In 
Sec. VI, some general remarks and applications to exponential product formulas 
are briefly mentioned. Summary and discussion are given in Sec. VII. 

II. Inner Derivation, its Formal Inverse and Uniqueness 

In the present section, we introduce the two hyperoperators (— 5~a~5b) and 
5 (-5-13)1 an d discuss the existence and uniqueness of these hyperoperators in the 
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domain T>a, which is defined by the set of convergent power series of the operator 
A in a Banach space. In general, it seems to be meaningless to use the symbol 
5 J 1 , because the inverse of the inner derivation 5 a does not necessarily exist and 
furthermore is not unique even if it exists. Fortunately in our problem, only the 
combinations (— 5^ 5 b) and 5^_ s -i B ^ appear in our quantum analysis of single- 
variable functions. Thus there is a possibility to define them uniquely. 

(i) Hyperoperator 5a-^b = —S^^b 

First we show that the hyperoperator (— 5 a 1 5b) is well defined when it operates 
on a function f(A) in the domain T>a- For this purpose, we confirm that 

Sa^bA = (-5 a 1 5b)A = 5 A \-5 B A) = 5 A l 5 A B = B, (2.1) 

namely 5a->b : A\—> B. More generally, we have 

n—l n—1 

5 A ^ B A n = £ A k {5 A ^ B A)A n - k - 1 = ]T A k BA n ~ k ~ l (2.2) 

fc=0 fc=0 

for any positive integer n. Thus, the hyperoperator 5a^b = —^^b is well 
defined, at least, in the domain T>a- Thus, the existence of 5a->b has been shown, 
but it is not unique. In fact, we put 

6 A ^ B f(A) = F(A,B) (2.3) 

which is constructed by the above procedure. Then, F(A, B) + G(A) may be also 
a solution of (—5 A 1 5 B )f(A), because 

- 5 B f (A) = 5 A F(A, B) + 5 A G(A) (2.4) 

for any operator G(A) in a Banach space. If we impose, besides the Leibniz rule, 
the linearity of the hyperoperator 5a-^Bi namely 

6a-+bU(A) + g(A)) = 5 A ^ B f(A) + 5 A ^ B g(A), (2.5) 
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and 

S A ^ B (af(A)) = a6 A ^ B f(A) (2.6) 

for a complex number a, then the uniqueness of 5a-^b is assured. In fact, the 
expression F(A, B) in (2.3) is obtained explicitly by using this linearity of the 
hyperoperator Sa-^b- 

In order to study the role of the hyperoperator 5 A -^b more explicitly, we intro- 
duce the symmetrized product {A m B n } sym (A ,b) by 

{A m B n } sym{A , B ) = E A kl BA k2 ■ ■ ■ A kn BA kn+1 , (2.7) 

fciH \-k n+ i=m,kj>0 

where m, n, {kj} denote non-negative integers. This symmetrized product is also 
written as 

(2.8) 



1 

{A m B n } sym{AjB) = 



-—(A + xB) m+n 



x=0 



Then, Eq.(2.2) is expressed by 

S A ^ B A m = {A rn ~ 1 B} sym (A,B) ■ (2.9) 

Hereafter, we write {• • -}sym(A,B) simply as {• • -} sym , when no confusion arises. 
Similarly we obtain 

S A ^ B {A m B} sym = {A m ~ l B% ym , (2.10) 

because 

m 

-S B {A m B} sym = -<fe(E A k BA m - k ) 

k=0 

m 

= E [(-S B A k )BA m - k + A k B{-5 B A m - k )} (2.11) 

k=0 

using the Leibniz rule. Using the commutativity of A and 5a and the relation 
(2.9), namely -5 B A m = 5 A {A m - l B} sym , we have 

m m—1 

-5 B {A m B} sym = Y,{$A{A k - 1 B} sym )BA m - k + E A k B(S A {A m - k - 1 B} sym ) 

k=l k=0 

= 5 A {A m - l B% yui . (2.12) 
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In general, we have the following formula. 

Formula 1 : For non- negative integers m(> 1) and n and for any operators A 
and B in a Banach space, we have 

- 5 B {A m B n } sym = 5 A {A m - l B n+l } sym (2.13) 

namely 

5 A ^ B {A m B n } sym = {A m ~ 1 B n+1 } sym (2.14) 

Consequently, the domain of the hyperoperator 5a->b is extended to the region 
^sym(A,B) which is a set of convergent noncommuting symmetrized power series of 
A and B. 

The proof of this formula is given as follows. First note that 

6a+xb(A + xB) m+n = 0, (2.15) 

namely 

- x5 B (A + xB) m+n = 5 A (A + xB) m+n . (2.16) 

By comparing the (n + l)-th terms of the both sides of (2.16) in x and using 
the relation (2.8), we obtain Eq.(2.13) and consequently Eq.(2.14). An alterna- 
tive derivation of Eq.(2.13) will be given by extending the procedure shown in 
Eqs.(2.11) and (2.12). 

Next we study the property of the power hyperoperators {S A ^ B }. It is easy to 
show the following formula. 

Formula 2 : For non-negative integers k, m(> k) and n and for any operators in 
a Banach space, we have 

6 k A ^ B {A m B n } sym = {A m - k B n+k } sym , (2.17a) 

and 

5 k A ^ B {A m B n } sym = if m<k. (2.17b) 



This gives the following general formula. 

Formula 3 : When f(A) is a convergent operator power series of an operator A 
in a Banach space, we have 

^ B {f(A)B m } sym = J o 1 dt 1 j\t 2 ---l^ 1 dtM {n \tnA)B m+n } 8ym 

= ^^J Q 1 dt(l-tr- l {f^\tA)B m+n } sym . (2.18) 

Here, f <yH \x) denotes the nth derivative of f(x). 

(ii) Hyperoperators S { _ s -i B) and d A ^ B = ^(s^b)-a 
An operator H defined by 

-5 A H = B (2.19) 

does not necessarily exist, as is well known. However, the hyperoperator 5^_ 5 -i B ^ is 
well defined, at least, when it operates on f(A) in the domain V A for an operator 
A in a Banach space. In fact, we can interpret it as 

m— 1 m— 1 

S { _ 6 - /B) A m = £ A k {5 { _ rA , B) A)A m - k - 1 = £ A k (5 A 5 A 1 B)A m - k - 1 

k=0 k=0 
m—1 

= £ A k BA m ~ k ~ l = 5 A ^ B A m . (2.20) 

k=0 

In other words, the formal hyperoperator 5 A X in S^_ s -i B ^ should be interpreted 
as a hyperoperator operating on the left-hand-side hyperoperator 5 A (not on the 
right-hand-side operator B). In this interpretation, the hyperoperator S^_ s -i B ^ is 
defined even when the operator (— S A X B) does not exist. 
In general, we obtain the following formula. 

Formula 4: Under the requirement of the linearity of the hyperoperators 5 a ^b 
and we have 

WW = 6 ^bI(A) = J l Q {^\tA)B} sym dt (2.21) 



for any operator f(A) G V A - Here, f <yl \x) denotes the first derivative of f(x). 

It should be remarked that the hyperoperator S^_ 5 -i dA ^ is a kind of differential 
defined only in the domain Va, whereas the hyperoperator 5a-*<ia is defined in 
a wider domain but is not a differential in the domain V sym (A,dA) outside of the 
domain Va- 

As was discussed before, the operator H = —S A X B does not necessarily exist, 
and it is difficult to define the power hyperoperators {^_§-i B ^} for n > 2 when 
H = —5~a~B does not exist. Furthermore, they are complicated 12 even if they 
do exist, unless H commutes with B. Thus, we define the following partial inner 
derivation 

d A ^B = ^_ s -i B y iA , {2.22) 

by which the commutator S^_ s -i B ^ is taken only with the operator A in a multi- 
variate operator f(A, B). For example, we have 

d A ^ B (ABA) = (8 { _ 6 - /B) A)BA + AB(S { _ 5 - /B) A). (2.23) 

This new hyperoperator (La-,b is defined in the domain T>a,b which is a set of 
convergent noncommuting power series of the operators A and B. Clearly, cIa^b 
is a kind of differential satisfying the Leibniz rule for B = dA. 

Next we study the power hyperoperators {d A _^ B }. Clearly they are also dif- 
ferentials defined in the domain V AB . It will be interesting to find the relation 
between d A _^ B and S A _> B . 

First note that 

(m— 1 \ 
E -WM"' ' 'J 

m-lj-l m-2m-j-2 

= E E A j ~ k ~ 1 BA k BA m ~ j ~ 1 + E E A j BA k BA m ~ j ~ k ~ 2 

3=1 k=0 j=0 k=0 

= 2{A m - 2 B% ym = 25\^ B A m (2.24) 
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for m > 2. In general, we obtain the following formula. 

Formula 5 : For m > n and for any operators A and B in a Banach space, we 
have 

d n A _^ B A m = n\5 n A ^ B A m . (2.25) 
We have also d A ^ B A m = for m < n. More generally, we have 

d n A ^ B f(A) = n\5^ B f(A), (2.26) 

when f(A) G V A . 

The proof of Formula 5 is given by mathematical induction using the following 
lemma and Formula 2. 

Lemma 1 : For non-negative integers m(> 1) and n and for any operators A and 
B in a Banach space, we have 

d A ^ B {A m B n } sym = (n + lMA™- 1 ^ 1 }^. (2.27) 

This is easily proved by using the definition (2.7) of {A m B n } sym as in Eq.(2.24). 
Formula 5 can be also confirmed directly from the consideration on the num- 
ber of permutations of B n . More intuitively, S A ^ dA denotes an ordered partial 
differential 11 . On the other hand, d\_^ dA denotes the nth differential, as will be 
discussed later. Consequently we have Formula 5. 

It should be remarked here that the hyperoperator d A _ >B is equivalent to 
8™_ s -i B ) when H = —d^B exists and it commutes with B. This equivalence 
has been already used implicitly in the previous papers 12 ' 13 . 

With these preparations, we discuss a general theory of derivatives of f(A) 
with respect to the operator A itself in the succeeding section. 
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III. Quantum Derivative, its Invariance and Operator Taylor Expan- 
sion 

In the present section, we give a general formulation of quantum derivatives 
{d n f (A) / dA n } which do not depend on the definition of the differential df(A). 
Our starting point of this general theory is that the differential hyperoperator u d v 
satisfies the Leibniz rule (1.9) and that it is a linear hyperoperator. 

(i) Quantum derivative and its invariance 
Now we start with the following identity 

Af(A) = f{A)A, (3.1) 

when f(A) 6 V A . Then, we have 

d(Af(A)) = d(f(A)A), (3.2) 

which is rewritten as 

(dA)f(A) + Adf(A) = (df(A))A + f(A)dA, (3.3) 
using the Leibniz rule. This is rearranged as follows : 

Adf(A) - (df(A))A = f(A)dA - (dA)f(A). (3.4) 

That is, we have 

5 A df(A) = 5 f{A) dA. (3.5) 

This is our desired formula on the differential df(A). 

In order to discuss the solution of Eq.(3.5), we rewrite Eq.(3.5) as 

5 A df(A) = -5 dA f(A). (3.6) 

Obviously, df(A) has a linearity property with respect to f(A). Thus, the solution 
df(A) of Eq.(3.6) is uniquely given in the form 

df(A) = 6 A ^ dA f(A), (3.7) 
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using the hyperoperator 5a—hIA = S A SdA introduced in Section II. This is also 
rewritten as 

df(A) = j l Q dt{f^{tA)dA] sym[AM) 

= £fM(A-t6 A )dt>dA, (3.8) 

using Formula 4, namely Eq.(2.21). 

The second equality of Eq.(3.8) is proven as follows. First we prove it when 
f(A) = A m for an arbitrary positive integer m. Clearly we have 

dA m = f o dt{fW(tA)dA} sym{AM) 

= (m j\ m - l dt){A m - l dA} syui[AM) 

= {A m ~ 1 dA} sym . (3.9) 

On the other hand, we obtain 

jT 1 f [l \A - tS A )dt ■ dA = jT 1 f {1 \(l - t)A + t(A - S A ))dt ■ dA 

in — 1 r i 

= m £ m ^C k (1 - t)H m - 1 - k dtA k (A - 5 A ) m - l - k ■ dA 
k=o u 

TO— 1 

= m E m -iC k B(k + 1, m - k)A k (A - 5 A ) m - l - k ■ dA 

k=0 

TO— 1 

= £ A k {A - 5 A ) m - 1 - k ■ dA 

TO— 1 

= £ A k (dA)A m ~ 1 ~ k = {A m ~ 1 dA} sym , (3.10) 

k=0 

using the beta function B(x, y), the binomial coefficient m Ck, the commutativity 
of A and 5 A , and the following relation 10 

(A-S A ) n -dA = (dA) • A n . (3.11) 

Thus, the second equality of Eq.(3.8) holds for f(A) G V A . 
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Furthermore we can derive the following relation. 

Lemma 2 : When f(A) G V A , we have 

6 f{A) = f(A)-f(A-6 A ). (3.12) 

Using this lemma, we obtain 

Sf(A) = f(A)-f(A-5 A ) 

= 6 A (fM(A) - -6 A fW(A) + ■■■+ { -^—8 n A 1 f in \A) + ■■■) 

= b A \j {X \A-tb A )dt. (3.13) 

This is formally written as 

H f®(A - tS A )dt = S?S m = S ~jf- ( 3 - 14 ) 

Thus, summarizing Eqs.(3.5),(3.7),(3.8) and (3.14), we obtain the following 
theorem on the differenial df{A). 

Theorem I : When f(A) G T> A , we have 

5 A df(A) = 5 f(A) dA, (3.15) 

and consequently 

df(A) = 5 A ^ dA f(A) 

= f 1 f {l \A-t5 A )dt-dA= 5 -^-dA (3.16) 
"'O o A 

for any choice of definitions of the differential df (A) . 

It should be noted that the ratio of the two hyperoperators 5f{A) an( A S A is well 

defined for f(A) G T> A , as was discussed in the preceding section. We define 10-12 

the quantum derivative df(A)/dA in Eq.(1.8), namely 

df(A) = d ?^.dA. (3.17) 
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That is, the derivative df(A)/dA is a hyperoperator which maps an arbitrary 
operator dA to the differential df(A) given by Eq.(3.16). Thus, we arrive at the 
following invariance theorem on the quantum derivative defined in Eq.(3.17). 

Theorem II (Invariance of Quantum Derivative) : When f(A) G Va, the 
quantum derivative df(A)/dA is invariant for any choice of definitions of the 
differential df(A) satisfying the Leibniz rule, and it is given by 

^tW„>^-W (3.18) 

Clearly, the ratio of the two hyperoperators Sf(A) an d Sa does not depend on 
the choice of definitions of the differential df{A). This invariance has been also 
discussed by Nazaikinskii et al 9 . in a different formulation based on the Feynman 
index method. The present confirmation of the invariance is more direct and 
transparent. 

(ii) Higher derivatives and operator Taylor expansion 

Now we discuss higher-order differentials {d n f(A)} and higher derivatives 
{d n f{A)/dA n }. 

(ii-1) Higher-order differentials and derivatives 

The hyperoperator dA-^B introduced in Eq.(2.22) is a derivation satisfying the 
Leibniz rule (1.9). Thus, dA-^B is a kind of differential hyperoperator, when 
B = dA. We prove here the following theorem. 

Theorem III: The nth differential d n f(A) is expressed by 

d n f(A) = d n A ^ dA f(A) (3.19) 

for any choice of definitions of the differential df(A), when f(A) G Va- 

The proof is given as follows. First note the following recursive formula 10 (3.21) 
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obtained by differentiating Eq.(3.5), namely 

5 A df(A) = 5 f[A) dA (3.20) 

repeatedly. 

Formula 6 : When f(A) G V A , we have 

6 A d n f(A) = n5 dn -i f{A) dA = -n5 dA d n - l f{A). (3.21) 

This gives the following result : 
Formula 7 : When f(A) G V A , we have 

d n f(A) = n5 A ^ dA d n - l f{A) = n\5 n A ^ dA f(A), (3.22) 

using the hyperoperator S a ^b = — 5~a~&b introduced in Sec. II. 

Here we have also used the relation df(A) = 5 A ^ dA f(A) given in Eq.(3.7). 
Using Formula 5, namely Eq.(2.26), we arrive at Theorem III. This result means 
that any differential hyperoperator d is generally expressed by 

d = d A ^ dA (3.23) 

in the domain V AjdA . 

Next we define 10 the higher derivatives {d n f(A)/dA n } through the relation 

d"f(A) = : dA ..... dA . (3.24) 

Here, d n f(A)/dA n denotes a hyperoperator which maps a set of the operators 

(dA, • • • , dA) =: dA -dA to d n f(A). In an ordinary mathematical notation, 

one may prefer to write as 

d n f(A) = ^V( dA,...,dA) . (3.25) 

n 
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However, as was emphasized before 10 , the product form (3.24) is essential in the 
present quantum analysis. That is, we use the product form (3.24) only when 
the derivative d n f(A)/dA n is expressed explicitly in terms of A and the inner 
derivations {5j} defined by 10 

8j : dA--- - - dA = {dAy-\5 A dA){dA) n - ] . (3.26) 

n 

If we use the notation (3.25), this property of product (3.26) and 

A: dA--- - ■dA = A{dA) n (3.27) 

n 

can not be shown explicitly. Clearly A and {5j} are commutable with each other. 

(ii-2) Integral representation of d n f(A)/dA n 

Here we express d n f(A)/dA n explicitly in an integral form in terms of the above 
inner derivations {Sj}. Our result is given by the following theorem. 

Theorem IV : When f(x) is analytic and f(A) G Va, any higher derivative 
d n f(A)/dA n exists uniquely for any choice of definitions of the differential df (A) , 
and it is given explicitly in the form 

^> = n\ j l d tl f dt 2 ... f;~ l dtnfM(A - ± tjSj). (3.28) 



.7=1 



Here, f^ n \x) denotes the nth derivative of f(x). 

The proof is given as follows. Once the above integral representation (3.28) is 
derived, the uniqueness of it is clear. In the case of n = 1, we have 

4?£) = *M = m - f} A - *A = Z 1 f^A - U A)d t (3.29) 
dA oa oa jo 

from Theorem II and Lemma 2. 

The nth derivative of f(A) divided by n\, namely f n (A, Si,---,5 n ) defined by 

1 d n f(A) 

f n (A, 6 h ---,6 n ): (dAT = ^ r j^ L : (dA) n = 5^ dA f(A) (3.30) 
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is shown from Formula 6 to satisfy the following relation 

(Si H h S n )f n (A, Si,--- , S n ) = f n -i{A, S h ■ ■ ■ , 6 n -i) - f n -i(A - Si, S 2 • • • , S n ). 

(3.31a) 

When f(A) = A k , Eq. (3.31a) means that 

5 A d n A k = n(d n - l A k dA - d n ~ l (A - 5 A ) k dA), (3.31b) 

which is equivalent by Formula 6 to saying that 

dAd n - 1 A k = d n ~\A - 5 A ) k dA. (3.31c) 

The solution of (3.31a) with the condition (3.29) for n = 1 is proven to be given 
by 

f n (A, £ n ) = /* dti Z' 1 dt 2 • • • jf*"" 1 <ft„/ (n) (^ - £ (3-32) 

i=i 

using the commutativity of ^4 and and the following formula for t = 1. 

Formula 8 : For any positive integers m and n, we have 

(xi + • • • + x n ) f dh f l dt 2 --- f n " dt n f( m+1 \tx - £ tjXj) 

i=i 

= / dti j dt 2 --- J dt n _i /<*») (t x - j2 t jXj ) 

^ i=i 

n— 1 \ 

J2t jXj+ i) , (3.33) 
i=i / 

when f(x) is a convergent power series of x. 

Proof of Formula 8 : Let h njTn (t; x,Xi,---, x n ) be the left-hand side of (3.33) 
minus the right-hand side of (3.33). Then, we have 

~J^h n ^ m (t, X, Xi, • • • , X^ = xh n ^ m +i(t, X, Xi, , X n ) + H n —i^ m (t, X Xi, X 2 , ' , X n ). 

(3.34) 

If we assume that /&n-i, m (t; x, xi, • • • , x n _i) = for all positive integers m and for 
any x, and {x 7 }, then we obtain 

~~^jh n ^ m (t , x, Xi, , x n ) = xh n ^ m -\-i(t, x, Xi, • • • , x n ). (3.35) 
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Thus, we derive 

d N 

^jyh"n,m(ti X, X\, , Xnj = xfl n ^ m +N{t, X, X\, ■ , X n ) (3.36) 

for any positive integer N. Thus, when f(x) is a polynomial of x, we have 
h n ,m+N{t', %, Xi, • • • , x n ) = for a large N. Clearly we have 

d k 



'h ni m{t, X, X\, , X n ) 



= (3.37) 

t=o 



dt k 

for any non-negative integer k(< N). The solution of Eq.(3.36) with (3.37) is 
given by 

h n , m (t;x,x h ■ ■ ■ ,x n ) = (3.38) 

for any positive integers n and m. Therefore, when f(x) is a convergent power 
series of x, we obtain Formula 8 by mathematical induction, because both sides 
of Eq.(3.33) is linear with respect to the function f(x). 

Thus, Theorem IV has been proven. An alternative proof of it is given in 
Appendix. The third proof is discussed in Sec. VI. 

(ii-3) Operator Taylor expansion and shift-hyperoperator Sa(B) 

Now we study the Taylor expansion of /(A + xB). First we prove the following 
general Taylor expansion formula. 

Theorem V : When f(A) G T>a, we have 

oo oo ry,n 

f(A + xB) = Y.^A^Bf{^)=Y.- ] d n A ^ B f{A) 

= V-":"7: ,; (3,9) 



n=0 n\ dA r 
Equivalently, 

f(A + xB) = S A (xB)f(A), (3.40) 
where the shift-hyperoperator Sa{B) is given by 

oo 1 

Sa(B) = £ -Ma-^bT = e dA - B . (3.41) 
n=0 n! 
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The proof of this theorem is given as follows. From Eqs.(2.8) and (2.17), we 
have 



1 



rn Am r m-ti nni 

°A^B /± — D Jsym — T 



d r 



dx 1 



■{A + xB) 



m 



x=0 



for m > n, and we have b n A _^ B A m = for m < n. Therefore, we obtain 



1 



d r 



dx 1 



-f{A + xB) 



\x=0 



for any positive integer n, when f(A) G T>a- This yields Theorem V. 
In particular, if we put B = dA, we obtain the following result. 

Theorem VI : When f(A) G T>a, we have 



OO r£ n 



f(A + xdA) = f(A) + £ ^d n f{A) = e xd f(A) 

n=l n\ 

with the differential hyperoperator d defined by (3.23) , namely 



d — d A - 



>dA- 



(3.42) 



(3.43) 



(3.44) 



(3.45) 



IV. Multivariate Quantum Analysis 

In this section, we formulate multivariate quantum analysis, in which we con- 
sider a set of noncommuting power series {f(Ai,- — ,A q )} = {f({Ak})}. This 
domain is denoted by T>{A k }^ namely f({Ak}) G ^{A k }- If we start from a complex 
number function /({#&}), it is a problem how to define the operator function 
/({Ah}), as is well known in quantum mechanics. Here, we start from the opera- 
tor function /({Afc}) itself which is specified in some appropriate procedures such 
as normal ordering. 

A definition of the partial differential djf({Aj~}) corresponding to Eq.(l.l) is 
given by 

d- f = lim /(Al> ' ' ' l Aj + MAj > '^AAJZ /( W) (41) 
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Norm convergence of Eq.(4.1) can be discussed in a Banach space and strong con- 
vergence is appropriate for unbounded operators. An algebraic partial differential 
corresponding to Eq.(1.2) is given by 

d j f({A k }) = [H j J({A k })] (4.2) 

with some auxiliary operators {Hj}. Both satisfy the Leibniz rule. In the present 
paper, we study general properties of multivariate quantum derivatives which are 
invariant for any choice of definitions of differentials. This invariance can be easily 
proved by extending the procedure shown in IIIA. Namely we have dj = Sb^a^Aj 
with Bj = — Sj\ The total differential df is defined by 

df=id j f = (£d j )f i (4.3) 

.7=1 3 

when / G ^{A k }- The nth differential d n f is also defined by 

d n f = (T,dj) n f. (4.4) 

3 

Clearly, {dj} commute with each other, namely djd k = d^dj, in the domain T*{A k }- 
One of the key points in the multivariate quantum analysis is to express d n f in 
the form 

d n f = n\ £ f\:l jn :dA h dA jn . (4.5) 

3^i'"i3n 

Then, we study how to calculate the hyperoperators {fj^...j n } in Eq.(4.5). 

(i) Ordered differential hyperoperator 

In order to study fj^...j n , we introduce here an ordered differential hyperoper- 
ator dj, ,•„... j as follows : 

dji,32,-,3n = (dj 1 dj 2 ' ' ' ^jn) ordered? (4-6) 

which means dj 1 j 2j ...j n f({A] c }) is given by those terms (found via the Leibniz rule) 
of dj x dj 2 • • • dj n f({Ak}) in which the differentials appear in the order dAj x dAj 2 • • • dAj r 
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For example, we consider an operator function f(A, B) = ABA 2 . Then we 
have 

d A , B f = (dA)(dB)A 2 , d B , A f = A(dB)(dA 2 ) = A(dB)[(dA)A + Ad A], 
d AA f = (dA)BdA 2 , d B , B f = 0. (4.7) 

Thus, using this ordered differential, we obtain the following formula. 
Formula 9 : In the domain T>{ Ak ^, we have 

dji---dj n =Y,d h ,..,j n . (4.8) 
p 

Here, Sp denotes the summation all over the permutations of (ji, • • • ,j n ). 

The proof will be self-evident. In particular, we have the following formulas. 

Formula 10 : In the domain T*{ Ak j, we have djdk = dj^ + dkj and 

d^ = n\d^. (4.9) 

n 

Formula 11 : In the domain T>{ Ak }, we have 

d n f = n\ £ d h ,., jn f (4.10) 

3li'"i3n 

for any positive integer n. 

Under these preparations, we find a procedure to calculate {/j"~.j n } in Eq.(4.5). 
In principle, they are obtained through the following relation : 

f\ n) , :dA n dA 1 =d il ... i f. (4.11) 

J 3li i3n 31 Jn Jl, ,JnJ V / 

Here, dj w ..j n f is expressed in the form 

dj u -,3nf = £ fkAdAjJhAdA^hj ■ ■ ■ (dA jn )f Kn (4.12) 
k 
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with some appropriate operators {fk,j}- In order to find {/)"... j n } explicitly, we 
have to rearrange Eq.(4.12) in the form of the left-hand side of Eq.(4.11). For this 
purpose, the following rearrangement formula 13 will be useful. 

Formula 12 (Rearrangement formula) : Any product Q1/1Q2/2 • • • Qnfn can 
be rearranged in the form 

n+l 

Ql/l • • • Qnfn = Yl /l/2 • • • fj-1 S dn{j,ji)dn{ji,j2) ' ' ' ^n(j k ,n+l) '■ Ql Qn (4-13) 

3=1 71 

with /o = 1. Here, denotes the summation all over the ways of the following 
division of the set of natural numbers (j, j + 1, • • • , n — 1, n) : 

(j,j + V--,™- = 7r(j,ji)7r(ji, j 2 ) • • -7r(jfc,n+ 1), (4.14) 

and 

7r(j, A;) = 0',i + V ••,£-!) (4-15) 
with j <ji < ••• < jk< n. Furthermore, the hyperoperator d^^) ls defined by 

using the partial inner derivation 5f.Q j = 8f-j which operates only on Qj in (4.13). 
The proof of Formula 12 is easily given by mathematical induction. 
It will be instructive to give here some examples : 

Qi/i = (/i " S fl ) : Q h 

Q1/1Q2/2 = (/l/2 - /lfy 2 ;2 ~ ^/i/ 2 ;l + fyi;lfy 2 ;2) : ' ^2, 

Q1/1Q2/2Q3/3 = (/1/2/3 - /l/2^/ 3 ;3 + /l^/ 2 ;2^/ 3 ;3 ~ /lfy 2 / 3 ;2 

~ S fifth;! + ^/i;i^/2/ 3 ;2 + ^/i/ 2 ;i^/ 3 ;3 ~ & frrf frrf f 3 fi) ■ Ql • Q2 ■ Qs- (4.17) 

(ii) Partial derivative and multivariate operator Taylor expansion 
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It will be convenient to define the following partial quantum derivative 

^ =n\&U, (4.18) 



OA, ■■■dA i . " 

using the hyperoperators j„} determined through the relation (4.11). Then, 
we obtain the following theorem. 

Theorem VII : When f({Ak}) G T>{ Ak j, we have 

oo 

/({Aj + XjdAj}) = E E x jl ---x jn d jl ... d J 

n=0ji,-,j n 

= E E x ii ' ' ' x jnfji,-,j n '• dAj 1 ■ ■ • dAj n 

OO ry . . . . ry . f)1l f 

= x * a :dA .... dAjAm 

n=0ji,-,j„ nl OA Jn 0n -Ji 

Equivalently, we have 

f({Aj + xjBj}) = exp ( E xjdA^n) f{{A 3 }). (4.20) 



In particular, 



^-3 ^3 

7=1 



OO rf n 



f({Aj + xdAj}) = e xd f{{A 3 }) = E ^d n f({A 3 }) (4.21) 

n=0 n! 



with d = Sj<ij. 



V. Auxiliary Operator Method 

It will be convenient to introduce the auxiliary operators {Hj} satisfying the 
following conditions: 

[H v H k ] = 0, [H j: A k ] = 0, [H h [H k , A k \] = 0, for j ^ k (5.1) 

and 

[Hj,[Hj,Aj}}=0. (5.2) 
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Using these auxiliary operators {Hj}, we introduce the following partial differen- 
tial 

djf^H^f^SnJ. (5.3) 

In particular, we have 

dA j = d j A j = [H j ,A j ], (5.4) 

and 

d 2 Aj = and dj(dA k ) = 0. (5.5) 
The total differential df is given by 

df = nHjJ] = (T,dj)f. (5.6) 

3 3 

One of the merits of this auxiliary operator method is that we can easily obtain 
the operator Taylor expansion as follows : 

3 

= f(i A j + XjShjAj}) = f({Aj + xjdAj}), (5.7) 
using Eqs. (5.2) and (5.4). That is, we have 

f({Aj + XjdAj}) = exp(E, x 3 S H j )f({A j }) = e E '^/({^}), (5.8) 

3 

using the relation (5.3). 

VI. Some General Remarks and Applications to Exponential Product 
Formulas 

It will be instructive to remark that when the operator A depends on a param- 
eter t, namely A = A(t), we have 10 

df(A(t)) = df(A(t)) dA(t) 

dt dA(t) dt { ' } 
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This formula insures again the invariance of the derivative df(A(t))/dA, because 
df(A(t))/dt and dA(t)/dt do not depend on the choice of the differential df(A(t)). 
Furthermore, we have 

df(g(A)) = df(g(A)) dg(A) 

dA dg{A) dA ' [ ' } 

because 

It will be also interesting to note the derivative of hyperoperators. The first 
differential of a hyperoperator /(Sa) is given by 10 

d(f(5 A )dA) = f Q dt[f^(t6i + 5 2 ) - /«(4 + tS 2 )} : (dA) 2 . (6.4) 

In general, we have 

d[f(A; 4, • • • , 5 n ) : (dA) n ] = £ dtfW(A - 6 2 ,---, 6 n+1 ) : (dA) n+l 

f( 1,k \A; 4, • • • , Sk-i,tkSk + 4+2, • • • , 4+i) 

{dA) n+l }. (6.5) 



k=l 



-f^ k \A; Si,---, 4-i, 4 + tkSk+i, 4+2, • • • , 4+i) 

Here, f (yl \x; Xi ■ • ■ , x n ) denotes the first derivative of f(x; Xi ■ ■ ■ , x n ) with respect 
to x and f^ 1,k \x; Xi ■ ■ ■ , Xk, • ■ ■ , x n ) denotes the first derivative of / with respect 
to Xk- Note that A and {4} commute with each other. These formulas will be 
also useful in proving Theorem IV. In fact, we obtain 

d 2 f(A) = d(df(A)) 

= d(j Q 1 dtf^\A-t5 A )dA) 

= J Q 1 dtd(f^(A-t5 A )-dA) 

= £ dh jf 1 dt 2 fW(A - tiSi - t 2 5 2 ) : (dA) 2 

+ jf 1 dt J Q \-s)ds[f^(A - s(tSi + 4)) - f {2 \A - s(Si + tS 2 ))} : (dA) 2 
= 2 jf 1 dh /J 1 dt 2 f^(A - tiSi - t 2 S 2 ) : (dA) 2 . (6.6) 
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Similarly we can derive Theorem IV using the above formula (6.5). 

There are many applications of quantum analysis to exponential product formulas 11-15 
such as the Baker-Campbell-Hausdorff formula. 

For example, if we put 

e Ai{x) e A i {x) mmme M') =e m j (6.7) 

the operator $(x) is shown to satisfy the operator equation 14 

^ = A-\*(x)) t exp(<5 Al(x) ) ■ • ■ exp((5. Vl(l) )A(^))^^ ( 6 . 8 ) 

using the quantum derivative of e : 

^ = °-^— = e * M - A) ., MA) = °— (6 . 9) 

The solution of Eq.(6.8) is given by 

log[exp(<^ l(f) ) • • • exp(^ r(t) )] 



j=1 'o exp(5 Al{t) ) ■ ■ ■ exp(5 Ar{t) ) - 1 



HA (t) 

x exp(tf Al(t) ) • • • exp( V l( t))A(^))— ^ + $(0). (6.10) 

This is a generalized BCH formula. 
In particular, we have 

log(eW) = / q _ ; log(e*VV^)J Aft + B. (6.11) 

Recursively we have 

log(e Al • • • e^) = fl dt ^ f ^ A, + E r (t)A r ) + $ 2 , r _ b (6.12) 

where $2,r-i = log(e A2 • • • e^- 1 ), and 

2? r (£) = exp(t6A 1 )exp(6A 2 ) • • • exp(^ r _ 1 )exp(^ j4r ). (6.13) 

The feature of these formulas is that $>(x) and log(e Al • • • e Ar ) are expressed 
only in terms of linear combinations of {Aj} and their commutators. 
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These formulas will be useful in studying higher-order decomposition formulas 



VII. Summary and Discussion 

In the present paper, we have unified an analytic formulation of quantum anal- 
ysis based on the Gateaux differential and an algebraic formulation of quantum 
analysis based on commutators, by introducing the two hyperoperators 6a-+b = 
—b~A&B and (Ia^b = ^(s^b^a- This general theory of quantum analysis gives 
a proof of the invariance of quantum derivatives for any choice of the definitions 
of differentials in the domain T>a- This domain can be easily extended 12 to the 
region T>a which is a set of convergent Laurent series of the operator A in a Ba- 
nach space. Multivariate quantum derivatives have also been formulated using 
the rearrangement formula. 

The present general formulation will be used effectively in studying quantum 
fluctuations in condensed matter physics and it will be also useful in mathematical 
physics. The present quantum analysis can also be extended to an infinite number 
of variables 14 . The quantum analysis has been also used 15 in extending Kubo's 
linear response theory 18 and Zubarev's theory of statistical operator 19 to more gen- 
eral nonlinear situations 11 . The invariant property of quantum derivatives derived 
in Sec. Ill is closely related 15 to the general feauture of the fluctuation-dissipation 
theorem 18-21 . General quantum correlation identities are also derived 15 using the 
quantum analysis. For the convergence of unbounded operators, see the second 
paper of Ref.17. 
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Appendix : Alternative Proof of Theorem IV 

First we study the case f(A) = A m for a positive integer m. The nth differential 
d n A m is expressed in the form 

d n A m = d^ dA A m = n\5^ dA A m = n\{A m - n (dAr} syui 

= n\ £ A ko (dA)A kl (dA)---A kn - 1 (dA)A kn 

kj>0,^2 kj=m—n 

= n\ £ A k °{A-5 l ) k ' 

... (A-5i 5 n ) kn :dA dA, (A.l) 

for m > n and d n A m = for n > m, using Theorem III, Formula 5, Formula 2, 
the definition of the symmetrized product, Eq.(2.7), and the following formula 10 . 

Formula A : For any operator Q, we have 

Qf(A) = f(A - 5 A )Q (A.2) 

when f(A) G V A . 

This yields Lemma 2. Now, we prove the following lemma. 

Lemma A : When f(A) = A m with a positive integer m(> n), the formula (3.28) 
holds. That is, we have 

E A k °(A - 5^ ■■■(A-5 1 5 n ) k - = F n (A; S 1 ■ ■ ■ , 6 n ), (A.3) 

kj>0,^2 kj=m—n 
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where 

771 1 r 1 /"ti rt — i ^ 

F n (A; *,...,(*„)= J / Q d*i / Q d* 2 • • • / Q n dt n (A - E *i^-) ro - n . (A.4) 



.7=1 



This lemma can be proved by mathematical induction as follows. We assume 
that Eq.(A.3) holds in the case of F n _i(A] 5i • • • , 8 n -i). Then, we have 

F n (A; 8 h ■ ■ ■ , S n ) = - m ' t dh /*' dt 2 --- /"*" 2 dt„_i 

(m — n + 1)! 7 o 70 jo 

1 



3=1 .7 = 1 

i- E <y B _ 2 )*»- 

fcj>0,^ kj=m— n+1 

x{(A -Si ^-i)"- 1 — (A — 6i 5 n ) k ^} 



5 n 



(A.5) 



under the assumption that Eq.(A.3) holds for F n _i(A; 5\, • • • , 6 n -i). Then, the 
above expression (A.5) can be rearranged as 

F n (A; 5 u ---,6 n ) = E A* (A - 5^ ■■■{A-5 l £ n _ 2 )*- 2 

kj>0,J2 kj=m—n+l 

x E {A -6! ^-i)^" 1 ■ (A - 8i 5 n f« 

<_i>o,A;>o,A!;_ 1 +fc;=A n _i-i 

E A k \A-5 1 f^--{A-5 l 5 n ) k ^ (A.6) 

fcj>0,^ kj—m—n 

by noting that fc + h + • • • + k n - 2 + + K = k + k\ + • • • + fc n _i - 1 = 
(m — n + 1) — 1 = m — n. Thus, we arrive at Lemma A. Any operator f(A) G V>a 
is expressed as a power series of {A m }. Then, Lemma A yields Theorem IV. 
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